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' Abstract. We report simultaneous observations in the X- 
I ray, infrared, and radio wavelengths of the galactic super- 
' luminal source GRS 1915+105. During episodes of rapid 
disappearance and follow up replenishment of the inner ac- 
cretion disk evidenced by the X-ray oscillating flux, we ob- 
. serve the ejection of relativistic plasma clouds in the form 
' of synchrotron flares at infrared and radio wavelengths. 
. The expelled clouds contain very energetic particles with 
' Lorentz factors of ^10^, or more. These ejections can be 
. viewed as small-scale analogs of the more massive ejecta 
' with relativistic bulk motions that have been previously 
! observed in GRS 1915+105. 

Key words: Accretion, accretion disks: Stars: individual: 
GRS 1915+105 - Stars: variables - Infrared: stars - Radio 
. continuum: stars - X-rays: stars 



1. Introduction 

The discovery of superluminal jets (Mirabel & Rodriguez 
1994) in the black hole X-ray transient source GRS 
1915+105 has opened the possibility of studying phenom- 
ena in our Galaxy that until recently were believed to 
be restricted to distant quasars and a few active galactic 
nuclei. In particular, it has been realized that since the 
characteristic dynamical times in the flow of matter onto 
a black hole are proportional to its mass, the events with 
intervals of minutes in a microquasar could correspond to 
analogous phenomena with duration of thousands of years 
in a quasar of 10^ Mq, which is much longer than a hu- 
man life-time. Therefore, the variations with minutes of 
duration observed in GRS 1915+105 in the radio, IR, and 



X rays sample phenomena that we have not been able to 
observe in quasars. 

X-ray observations of GRS 1915+105 with the Rossi X- 
Ray Timing Explorer (RXTE) revealed remarkable quasi- 
periodic oscillations (QPOs) that are believed to arise in 
the accretion disk around a black hole of stellar mass 
(Greiner et al. 1996, Morgan et al. 1997; Belloni et al. 
1997a,b; Chen et al. 1997). Among the diversity of QPOs 
observed in GRS 1915+105 there is a class of recurrent 
episodes with amplitude variations of ~10^^ erg s^^ (at 
a distance of 12.5 kpc) in the X-rays, in time scales of 
one minute to tens of minutes. These QPOs have been 
attributed to the rapid disappearance and slower replen- 
ishment of the inner region of the accretion disk (Belloni 
et al. 1997a, b). If the accretion is an advection-dominated- 
fiow (Narayan et al. 1997), the disappearing mass will go 
quietly though the horizon of the black hole. Furthermore, 
if a fraction of the mass of the accretion disk were blown 
away, one should see synchrotron emission from expanding 
clouds in the radio, and perhaps shorter wavelengths. 

Although previous observations have shown flares with 
similar periodicities in the radio (Pooley & Fender, 1997; 
Rodriguez & Mirabel, 1997a) and infrared (Fender et al. 
1997), no truly simultaneous observations in the X-ray, in- 
frared, and radio wavelengths have been reported so far. 
The observations reported here firmly establish the gene- 
sis of expanding clouds of relativistic plasma when GRS 
1915+105 recovers from large amplitude dips in the X-ray 
flux. 

2. Multiwavelength observations and light curves 

In this Letter we report observations made on May 15 and 
Sept 9, 1997. On May 15 at the UT intervals 6-10 h and 
12-16 h, we used the 27 antennae of the Very Large Array 
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(VLA) split in three sub-arrays to allow simultaneous ob- 
servations at 2, 3.6, and 6 cm. On Sep 9,17 antennae of the 
VLA were available for this project, and the observations 
were done at 3.6 cm only. 

The infrared observations were carried out with the 
United Kingdom Infrared Telescope (UKIRT). On each 
date, sixty images with the IRCAM3 256x256 InSb detec- 
tor array, with a K-band filter of 0.37/xm width centered 
at 2.2/xm, were obtained continuously every 1 min for an 
interval of h. The 2.2/im flux of GRS 1915 + 105 
on time bins of 1 min was computed using relative pho- 
tometry with several stars in the field and calibrated with 
well known standard stars. The infrared fluxes were dered- 
dened assuming Ax =3. 3 mag (Chaty et al. 1996). 

The X-ray observations were done with the PGA on 
RXTE. It consists of five nearly identical proportional 
counter units with detectors sensitive to X-rays between 
2 and 60 keV. 

The results for May 15 arc shown in Figure 1. The 
RXTE observations exhibit large amplitude variations 
with sudden drops of the flux into dips of ~5000 counts 
that last for at least 40 min. The recovery of the X- 
ray flux from the dips sometimes ends in sharp peaks that 
reach ~10~^ erg cm~^ s~^ (2 mJy) with typical durations 
of ^25 s. These rebounds of the X-ray flux are followed 
by smaller amplitude oscillations of ~60 s period that last 
tens of minutes. At the time of an X-ray gap around 14 h 
UT, an infrared flare of 12 mJy was observed. The rise of 
the X-ray flux observed after 14.25 h UT suggests a dip 
during that gap. Interpolating the pattern, a sharp peak 
and sequences of oscilations would have occured at ~13.7 
h UT. 

The quasi-periodic oscillations in the radio flux are 

typical of the optically thick radio core state, usually ob- 
served when the source oscillates between 10 and 100 mJy. 
A time-shift of the flux peaks, with the short wavelengths 
peaking flrst, is clearly observed. 

The results for Sep 9 are shown in Figure 2. The PGA 
count rates show large amplitude oscillations (~10^® erg 
s""'^) every '-~^50 s for intervals of tens of minutes, which 
are followed by sudden drops into dips. The recovery from 
the dips is slower than the time of decay and lasts a few 
tens of minutes. In the X-ray light curves of May 15 and 
Sep 9 we note the following similarities: 1) both the ~60 s 
oscillations on May 15 and the ~50 s oscillations of Sep 9 
disappear during the dips, and 2) the low intensity vahie 
in both cases is about 5000 counts s~^, which is well above 
the background rate of 100 counts s~^. However, there are 
obvious differences between the light curves in the two 
dates. 1) On May 15 the X-ray dips lasted longer than 
on Sep 9. 2) On May 15 the infrared flare had smaller 
amplitude and duration than the radio flare, whereas on 
Sep 9 they had similar amplitude and duration. 3) On 
May 15 the time-shift between the infrared peak and the 
following peak at 3.6 cm is 40 min, on Sep 9 it is 16 min. 



The relation between the X-ray, infrared, and radio 

light curves for Sep 9 is shown in more detail in Figure 3. 
The X-ray spectral index (13-60 kev)/(2-13 keV) suddenly 
increased at the beginning of the dip. Despite the 10^ or- 
der of magnitude difference in wavelength, the infrared 
flare and the follow up radio flare are strikingly similar 
in amplitude, duration, and slope of the increasing and 
decreasing fluxes. The infrared flare seems to start during 
the recovery from the X-ray dip at ~8.23 h UT, when the 
X-ray spectrum softened again and an X-ray spike was 
observed. The rise of the infrared flux to maximum value 
continues for a few minutes after the 50 s oscillations in 
the X-ray flux appeared again. 

3. Infrared synchrotron precursors 

The radio emission from GRS 1915+105 during 1997 May 
15 shows a behavior consistent with that expected for syn- 
chrotron emission from an adiabatically expanding cloud 
(van der Laan 1966). In all four bursts observed in the 
period of 6 to 16 h UT, the light curves at 6, 3.6, and 2- 
cm reach their peaks with the characteristic wavelength- 
dependent delay (the short wavelengths peaking first). 
These delays rule out the possibility that the observed 
events were produced by time-variable free-free opacity, 
as proposed by Rodriguez & Mirabel (1997a) to explain 
sinusoidal variations observed in 1995, since in this model 
no wavelength-dependent delays are expected. When the 
formation of clouds is faster than the decay of the flux, 
consecutive radio events blend together and the integrated 
emission from GRS 1915+105 may appear with a sinu- 
soidal modulation (Rodriguez & Mirabel, 1997a). 

In the van der Laan (1966) model, we can estimate 
p, the energy spectral index of the relativistic electrons 
{N{E) oc E~P), from the equation that relates the ob- 
served maximum flux density, Sm,\, at a given wavelength, 
i.e, Sm,\ oc X-(.'^P+^)/('^P+'^). For this estimate we use the 
last burst observed on May 15 between 14 and 16 h UT, 
since it appears to be more isolated than the others. Since 
Smficm = 39 mJy, and Sm,3.6cm = 51 mJy, we obtain p ~ 
0. 

The ejection is defined to occur at t = to- The time 
tm,\ since ejection in which the light curve at a given wave- 
length reaches maximum is given by tm.x oc \iP+^) / i^P+^) . 
With p ~ we obtain tQcm/tsscm = 1-4. Furthermore, 
from the data in Figure 1, we have iecm + ^o = 14.95 h UT 
and tsscrn + to = 14.70 h UT. Therefore, teem - tsscm = 
0.25 h, and iecm — 0.9 h. We then find that the ejection 
of the plasma occurred at to = 14.05 h UT. 

In this simple model, the UT time of maximum 
flux density at a given wavelength, for the particular 
event analyzed above, is given by tm,\{UT h) = 14.05 + 
0.9(A/6cm)(2/3). Then, the maximum flux density at short 
wavelengths (i.e. the near infrared) are observed very 
shortly after the ejection. It is only in the radio wave- 
lengths that significant time delays occur. We note that 
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Fig. 1. Light curves of GRS 1915 + 105 for May 15, 1997 obtained with the VLA at 2, 3.6, and 6 cm, with UKIRT 
at 2.2/im, and the PCA of RXTE at 2-60 keV. The time-lag of the peaks as a function of wavelength is the time 
scale for the clouds to expand and become transparent at longer wavelengths. It is shown that the infrared peak at 
~14 h UT occurred at the time expected for a synchrotron precursor of the peaks at radio wavelengths appearing 
later in the interval 14.6-14.9 h UT. The infrared flux densities have been dereddened by Ax=3.3 mag. The units 
on the left of the vertical axis correspond to the flux densities at radio and infrared wavelengths, those on the right 
to the X-ray counting rate. Typical error bars for each wavelength are shown at the top right corner. The gaps in 
the X-ray flux every ~90 minutes are due to Earth occupations and South Atlantic Anomaly passages. 
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Fig. 2. Light curves of GRS 1915 + 105 for Sep 9, 1997 
obtained with the VLA at 3.6 cm, with UKIRT at 2.2/im, 
and the PCA of RXTE at 2-60 keV. A full oscillation of 
the infrared synchrotron precursor of the following radio 
outburst is observed. Note the striking similarity of the 
flares at infrared and radio wavelengths. The units on the 
left of the vertical axis are for the radio and infrared flux 
densities, whereas those on the right are for the X-ray 
counting rate. Typical error bars for each wavelength are 
shown at the top right corner. 
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Fig. 3. Blow up of the light curves at the time of the 
infrared flare detected on Sep 9, 1997. The infrared flare 
starts during the recovery from the X-ray dip, when an 
X-ray spike was observed. These observations show the 
connection between the rapid disappearance and follow 
up replenishment of the inner accretion disk seen in the 
X-rays, and the ejection of relativistic plasma clouds 
observed as synchrotron emission at infrared and radio 
wavelengths. The hardness ratio (13-60 keV)/(2-13 keV) 
is shown at the bottom. 
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the 2.2/um emission indeed reached maximum near 14.0 

h UT. This result strongly supports the interpretation of 
the IR burst on May 15 as the synchrotron precursor of 
the radio bursts. To our knowledge, this is the first time 
in X-ray binaries that a clear association is established 
between an IR synchrotron precursor and the follow-up 
radio bursts. 

4. Jet formation during the X-ray dips 

Figure 3 shows that on Sep 9 the infrared flare started 
during the rise of the X-ray flux from the dip, when a 
softenning of the spectrum and an X-ray peak were ob- 
served. Since at infrared wavelengths the ejected plasma 
becomes transparent very shortly after the ejection, the 
rise of the infrared flux to maximum indicates that the 
injection of relativistic particles lasted '--^10 min, or less. 

If we assume that the synchrotron emitting plasma 
expands at ~0.2c, as observed in a larger scale ejection 
(Mirabel & Rodriguez 1994), 15 minutes after the ejec- 
tion the clouds have dimensions of ~ 10^"^ cm. Using this 
dimension and the infrared and radio flux density values 
measured on Sep 9, with equipartition arguments we es- 
timate a brightness temperature of 10^^ K, a magnetic 
field of 16 G, an equipartition energy in relativistic elec- 
trons of 5 X 10^^ erg, and a typical synchrotron luminosity 
integrated from 2.2fj,m to 6 cm of the order of 10'^^ erg 
s~^. These parameters are comparable to those derived by 
Fender et al. (1997) from observations on other epochs. If 
the electrons that produce the synchrotron radiation have 
a representative Lorentz factor of ~10^, and there is one 
proton per electron, the minimum mass of the clouds that 
are ejected every few tens of minutes is ~ 10^^ g. 

For the parameters of these clouds the energy losses are 
dominated by adiabatic expansion and the clouds should 
have a lifetime of about 1 h, as observed here. An in- 
frared jet as the one reported by Sams, Eckart & Sunyaev 

(1996) , which reached distances of ~0.2 arcsec (2500 AU 
at 12.5 kpc) requires very energetic particles, and a re- 
acceleration mechanism. We point out that radio blobs 
with bulk motions of 0.92c and lifetimes of several weeks 
have been observed (Rodriguez & Mirabel, 1997b). 

5. Conclusions 

1) At times of quasi-periodic oscillations of large ampli- 
tude in the X-rays, relativistic plasma clouds emanate 
from the X-ray binary. The time shift of the emission at 
radio wavelengths is consistent with synchrotron radiation 

from gas in adiabatic expansion. 

2) The association of infrared synchrotron precursors to 
follow-up bursts at radio wavelengths is firmly established. 

3) By the simultaneous observations reported here, we 
definitively confirm that the infrared and radio flares are 
associated to X-ray dips, as suggested by Fender et al 

(1997) and Pooley & Fender (1997). The onset of the ejec- 



tion takes place during the recovery from the X-ray dip of 

hard (13-60 keV)/(2-13 keV) spectrum, at the time of a 
spike in the X-rays. 

4) The expanding clouds of plasma reported here can be 

viewed as small-scale analogs of the more massive ejections 
associated to large radio outbursts (Rodriguez & Mirabel, 
1997b). The equipartition energy in relativistic electrons 
in the mini-ejecta is of order of 5 x 10^^ erg, namely, ^10^ 
times less energetic than in major ejecta. However, the 
genesis of these smaller scale clouds is much more frequent, 
and may represent an important phenomenon in the long 
term evolution of the binary and its environment. 

5) The multiwavelength variations with minutes of dura- 
tion in GRS 1915-1-105 show how the sudden disappear- 
ance of the inner accretion disk around a black hole of 
stellar mass triggers the formation of relativistic expand- 
ing clouds. Analogous phenomena in quasars with black 
holes of 10^ M0 would require thousands of years of obser- 
vations. Therefore, microquasars in our own Galaxy can 
provide new information on the physics of accretion disks 
around black holes and are thus useful to gain insight into 
the relativistic ejections seen elsewhere in the Universe. 
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